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We study the capabilities of the Majorana DEMONSTRATOR, a neutrinoless double-beta decay experiment 
currently under construction at the Sanford Underground Laboratory, as a light WIMP detector For a cross 
section near the current experimental bound, the MAJORANA DEMONSTRATOR should collect hundreds or even 
thousands of recoil events. This opens up the possibility of simultaneously determining the physical properties 
of the dark matter and its local velocity distribution, directly from the data. We analyze this possibility and find 
that allowing the dark matter velocity distribution to float considerably worsens the WIMP mass determination. 
This result is traced to a previously unexplored degeneracy between the WIMP mass and the velocity dispersion. 
We simulate spectra using both isothermal and Via Lactea II velocity distributions and comment on the possible 
impact of streams. We conclude that knowledge of the dark matter velocity distribution will greatly facilitate 
the mass and cross section determination for a light WIMP. 

PACS numbers: 95.35.+d 



I. INTRODUCTION 

Direct detection experiments offer the possibility of a non- 
gravitational detection of the most common form of matter in 
the Universe 1 1 1. Swift progress is being made in this field, for 
example the LUX and XENON IT experiments are expected 
to push the bound on the spin-independent cross section of a 
100 GeV WIMP to the unprecedented 10""'^ cm^^ level. Per- 
haps an even more dramatic development, however, may come 
about on the fight WIMP front (< 10 GeV). There, the im- 
provement relative to the current bound may be even greater, 
as a result of the MAJORANA DEMONSTRATOR lE) 13, as 
shown in Fig. [T] below. If the light WIMP cross section lies 
near its current bound, the MAJORANA DEMONSTRATOR is 
poised to collect hundreds or even thousands of recoil events, 
opening up an interesting possibility of not only discover- 
ing the dark matter particle, but also accurately measuring its 
properties. We will explore this possibility, and the issues that 
arise in connection with it, in this paper 

Before setting up our analysis, let us note that such figfit 
WIMPs are not without theoretical or experimental motiva- 
tion. On the theoretical side, 1-10 GeV dark matter (DM) is 
typical in models relating the baryon and DM abundances |4- 
[TOJ (see m il and references therein for recent work). In these 
models, the DM has a nonzero particle-antipaiticle asymme- 
try, thereby suppressing the indirect signatures of DM anni- 
hilation (though see |12 | for a counterexample). Establisfiing 
sucfi low-mass DM at a collider will be difficult since although 
the Tevatron and the LHC can detect the missing energy asso- 
ciated with light DM, they lose sensitivity to the mass of DM 
when it is much less than 100 GeV (see e.g., 1 13 1). The miss- 
ing energy could then be attributed to a variety of sources, 
from novel neutrino interactions [14| to extra-dimensional 
particles lfTSl [T6|. Hence the fate of this well-motivated class 
of models may hinge on a low-threshold direct detection ex- 
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periment such as the MAJORANA DEMONSTRATOR. 

On the experimental side, the interest of the community has 
been piqued over the last years by the CoGeNT and DAMA 
experiments, which have claimed signs of light dark matter 
ifTTl [TSl . With its remarkable sensitivity, the MAJORANA 
Demonstrator is expected to resoundingly refute or con- 
firm this possibility. In the latter case, it will see tens of thou- 
sands of events. In the latter case, it will be in agreement with 
the results of the null experiments |fT9]j22l . 

Let us suppose dark matter is discovered at the MAJO- 
RANA Demonstrator. How well could one infer its mass 
and cross section and what obstacles would one encounter? 
Thinking of this as a problem with a beam and a detector, one 
obvious answer is: the uncertainty of the beam. "The beam" in 
this case is provided by the local dark matter distribution and 
its characteristics are not uniquely set, as we discuss below. To 
avoid any confusion, we stress that for the present purpose the 
astrophysics obviously cannot be "integrated out", along the 
lines of ||23] |24l (see ll25l and ll26l for applications to recent 
experiments). The power of the "integrating out" technique is 
that it allows mapping the results of one experiment into an- 
other in an astrophysics-free manner. To determine the parti- 
cle physics properties of dark matter from a given experiment, 
however, the astrophysics needs to be "integrated in". 

Perhaps the answer lies in the large number of events: 
Given large statistics, can one directly measure both the fun- 
damental properties of DM as well as its local phase space dis- 
tribution? We will aim to provide a quantitative answer to this 
question, under "the most vanilla" of assumptions: that the 
DM has an elastic spin- and momentum-independent interac- 
tion, and that the DM is distributed in the Maxwell-Boltzmann 
fashion. We make these assumptions deliberately, as they 
form the standard framework any putative signal is likely to be 
analyzed under, at least at first. More precisely, we generate 
our mock signals with both isothermal and Via Lactea II nu- 
merical distributions (we also comment on the possible impact 
of the stream substructures). It is the fit analysis of the mock 
data that we deliberately constrain to the "the most vanilla" 
framework. Furthermore, we will remain agnostic about pre- 
viously mentioned putative signals of DM detection lilTl ITSll 
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(which would yield even larger signals), and will assume that 
the DM may be just at the border of what is allowed by the null 
results (see the black curves in Fig.[T]for current constraints). 

We draw attention of the reader to several important earlier 
studies that have examined the effect of astrophysical uncer- 
tainty on mass and cross section determination Il27l - l32ll . Here 
we extend these results to low masses and specifically focus 
on the capabiHties of the Majorana Demonstrator. We 
also outline an important degeneracy that arises when one si- 
multaneously fits the astrophysics and the dark matter proper- 
ties and show that it impacts the determination of the WIMP 
mass in an essential way. 

In the next section we review the basic theory of direct de- 



tection. In Sect. Ill we discuss our fit methodology and mock- 
up of the Majorana Demonstrator. In the first subsec- 



tion of Sect. IV we present the results of fits to the mass, cross 
section and the local velocity dispersion, assuming the "true" 



(input) velocity distribution is isothermal. In subsection IV B 



we examine to what extent light DM is sensitive to differ- 
ent forms of the velocity distribution by inputing the distribu- 
tion from the high-resolution N-body simulation. Via Lactea 
II. In Sect. IV C we comment on the degeneracy at heavier 
masses (m^ ^ 100 GeV) and in Sect. IV D we briefly touch 



on the signatures of streams. We summarize our conclusions 
in Sect.E 



II. BASICS 

While the nature of WIMP-nucleus scattering is obviously 
unknown, in the vast majority of existing scenarios the scat- 
tering is mediated by a sufficiently heavy mediator particle 
(compared to the momentum transfer in the direct detection 
process). Upon integrating out the exchanged particle, in 
the usual effective field theory approach, one finds that all 
the physics of the DM-nucleus scattering is contained in a 
set of higher dimensional operators. Several studies of the 
range of possibilities have been earned out |33-35|. A com- 
plete analysis of DM direct detection would in principle ex- 
amine how each of these operators is sensitive to astrophysi- 
cal uncertainty in turn - a daunting, but eventually necessary 
task. In the interests of clarity of the presentation, we will 
vary the astrophysics while restricting ourselves to the sim- 
plest and most-studied interaction form: isospin-conserving, 
spin-independent scattering that depend on neither the incom- 
ing DM velocity nor the exchanged momentum. Such an cross 
section could come from a scalar interaction {qq){XX) (such 
as from Higgs-exchange for a neutralino, e.g. ||36l ) or a vec- 
tor interaction {li^i^i_q){X^^^X) (arising from the exchange of 
a Z' vector boson). Our approach is similar to that of l,28„ 29j 
in the case of future signals and ll37l in the case of exclu- 
sions from current experiments. A complementary study was 
carried out in |38| where there the authors instead fixed the 
astrophysics and considered a discrete set of particle physics 
choices: with different types of momentum-dependent scat- 
tering and the standard momentum-independent scattering. 

To obtain the average differential rate per unit detector mass 
of a WIMP of mass mx scattering on a target nucleus of mass 



rriN, one convolves the cross section with the DM velocity 
distribution (see Ii36ii39i for reviews) 
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where Ve {t) is the velocity of the laboratory observer with re- 
spect to the galactic rest frame, f{v) is the DM velocity dis- 
tribution in the rest frame of the galaxy, and for elastic scat- 
tering, Vmin{EFi) = \/mAri?fl/2/i^ is the minimum DM ve- 
locity in the lab frame to produce a nuclear recoil of energy 
En. We have introduced the DM-nucleus reduced mass jix- 
We adopt the standard fiducial value of the local DM density 
Po = 0.3 GeV/cm^. For spin- and momentum-independent 
interactions, the differential cross section can be written as 



da 



dEn 2^^%v^ 



-^TZTi^SI^ \Er,A), 



(2) 



where the nuclear form factor F{Eji^ A) accounts for the mo- 
mentum transfer to the nucleus. The Helm form factor we use 
can be found in ll40l . The spin-independent DM-nucleus cross 
section is 
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where /i„ is the DM-nucleon reduced mass, and fp and /„ 
are the DM couplings to protons and neutrons respectively. 
Throughout, we will make the standard simplifying assump- 
tion of isospin-conserving scattering, /p = /n = 1. With 
these assumptions, the scattering rate simplifies to 
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where we have introduced the mean inverse speed 
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With this framework the unknown quantities are the DM mass 
mx, scattering cross section on a nucleon (7„ and the local 
velocity distribution f{v) or equivalently g(w„im))- We now 
discuss possible forms of the velocity distribution. 



A. Velocity Distributions 

In the first estimates of the rate of DM-nucleus scattering, 
Goodman and Witten simply used the average velocity of a 
DM particle in the halo Since then, a canonical frame- 
work has emerged in which the DM halo is taken to be an 
isothermal, Maxwell-Boltzmann (MB) distribution ||4T] l42l . 
with a cut-off corresponding to the escape speed. In the galac- 
tic rest-frame the distribution. 
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V < Vesc 
V > Vesc, 
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is fully specified by just two parameters: its dispersion vq and 
the local escape speed v^sc- The relevant velocity integral 
Eq.(j5]) and the normalization constant N have a closed form 
expression (see Appendix). In the first stage of discovery a 
signal is likely to be fit under this simplest and most studied 
assumption, which is why we choose it for our analysis. 

Throughout, we will examine time-averaged values, adopt- 
ing (fe(i)) = 230 km/s |43, 44]. The dispersion, vq, is less 
well known. In principle, given a model of the Galaxy, it could 
be related to observable properties of the DM halo, as done in 
for example EtII . To see how well the experiment can on its 
own determine the value of vq, we leave it as a free parameter 
In some calculations below, we do impose priors of various 
strength to gauge how the uncertainty of vq impacts the deter- 
mination of (T„ and mx ■ For the local escape speed we adopt 
the fiducial value 544 km/s, but note that the recent RAVE 
survey finds that at 90% CL the escape speed is in the region 
498 km/s < Vesc < 608 km/s Eg. 

Detailed N-body simulations appear to be better fit by NFW 
profiles, the velocities of which show non-Maxwellian fea- 
tures, especially on the high end (46]. We will investigate 
in Sect. |IVB| if Majorana Demonstrator, given mas- 
sive statistics, would be able to tell a difference between 
the (averaged) numerical velocities and the isothermal model. 
Even more distinctive for direct detection is the possibil- 
ity of unequilibrated velocity structures coming from the re- 
cent tidal destruction of DM subhalos. These can include 
streams fl7ll48l . as well as more spatially homogeneous tidal 
debris in the form of sheets and plumes P9l [5P|. We will 
discuss some implications from streams in Subsec. [rVD| 



III. EXPERIMENTAL DETAILS AND FIT SETUP 

Though the MAJORANA DEMONSTRATOR f3\ is primarily 
a neutrinoless double-beta decay experiment, the experimen- 
tal setup is well-suited for DM direct detection as well fj). 
The Demonstrator plans to achieve a very low threshold, 
in electron-equivalent units, 0.3-0.5 keVoc [2.1 (or equivalently 
Enji « 1.4 — 2.2 keV in nuclear recoil). Such low thresholds 
are due to P-type contact detector technology f5J|, also used 
by the CoGeNT experiment. Here we will take 0.3 keVoo to 
compare with [2]. The collaboration plans to have roughly 
40 kg of Germanium detectors, deployed in three stages over 
the next two years. Here we will assume an exposure of 100 
kg-yr [2|. The DEMONSTRATOR is now being consti'ucted at 
the Sanford Underground Research Facility (SURF) in Lead, 
South Dakota |[3l- They plan to being data-taking in 2013, 
with the first set of two strings. 

We also model the Xenon experiments, LUX and 
XENON IT in a simplified manner LUX is currently under 
construction at SURF while XENONIT will be installed at 
Gran Sasso National Laboratory in Italy, with planned data 
taking to begin in 2015. We mock up the XENONIT ex- 
periment by assuming a 2.2 ton-yr exposure with an energy 
window 8.4 — 44.6 keV and a 40% acceptance ll54l . We have 
checked that for heavy masses this mock-up agrees reasonably 
well with Fig. 2 of |55|. For the LUX experiment, we take 
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FIG. 1: Sensitivity reach of the upcoming and current experi- 
ments. The current constraints are shown with black curves: 2012 
XENON 100 |21| along with the low-threshold analyses done by 
XENONIO |20|, CRESSTl |51| and DAMIC |52|. For the upcom- 
ing experiments, we use colored curves: solid assuming the veloc- 
ity dispersion vq — 230 km/s and dashed for vo ~ 170 km/s and 
vo = 290 km/s. For the MAJORANA DEMONSTRATOR we assumed 
a 0.3 keVee energy threshold and 100 kg-yr exposure. For LUX, we 
took a 100 kg-yr exposure with Eth ~ 5 keV, and for the XENONIT, 
a 2.2 ton-yr exposure and Eth =8.4 keV (see |IV Cj for more details). 
The black circles indicate points this paper examines in detail. 



30000 kg-day exposure, the recoil energy window 5-30 keV 
and 45% efficiency. In both Xenon experiments, our analysis 
is deliberately conservative in that we calculate signal spec- 
tra in a fixed energy window, and do not convert nuclear re- 
coil energy into of expected number of scintillation photons 
and ionization electrons. Doing so effectively lowers the en- 
ergy threshold of the Xenon based experiments, but to do so 
properly requires a more detailed understanding of the exper- 
iments than we are capable of yet. Thus in comparing with 
the XENONIT limits in ll55l for example, our results repro- 
duce well the exclusion sensitivity at high mass but are less 
stringent than what appears in [55 1. We urge the Xenon col- 
laborations to improve on the analysis presented here. 

The DM sensitivity of MAJORANA DEMONSTRATOR is 
crucially affected by the background from cosmogenically ac- 
tivated tritium. This has been well studied in the thesis |2|. 
Assuming that a detector spends a total of 15 days at the sur- 
face, the average background rate contributed by decaying 
tritium is 0.03 counts/day /kg/keV |2|. The tritium spectrum 
is well-known and extends to around 18 keV.e- This back- 
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ground can of course fluctuate up and fake a WIMP signal. 
We will assume that this is the only relevant background in the 
signal region, and note that the presence of unknown "surface 
events" would significantly comprise their sensitivity (see for 
example |56|). 

For our statistical analysis, we bin events in Eji (taking 0.9 
keV bins for Majorana Demonstrator, and 3 keV bins 
for the Xenon experiments). The standard prescription for the 
likelihood of parameters is a product of Poisson distributions 
across the energy bins: 



m = n 



(7) 



where Nobs{Ei) is the observed number of events in the i*^ 
energy bin and N{6,Ei) is the number of expected events 
from parameters 9 in the same bin. Given that at MAJORANA 
DEMONSTRATORwe will be examining cases with large num- 
bers of events (from the signal plus the background), the Pois- 
son distributions actually become Gaussians and the overall fit 
becomes well-described by a standard statistics. We have 
verified this explicitly. 

Moreover, since we are interested in the results obtained on 
average at a given experiment, we follow the method outlined 
in the Appendix of |57|, and calculate the average 



(8) 



where bi and Si are the background and signal expectation in 
the i*'' energy bin. In the case of an exclusion, we calculate the 
limit by imposing (x^) < acLNdof, where acL depends on 
the desired confidence level. Throughout we will use 95% CL. 
When producing confidence regions from input signal spectra 
we calculate £ cx e^^-*^ at each parameter point, and then 
integrate the likelihood over the subspace of parameters that 
yield 95% of the total likelihood. This method is useful for the 
non-Gaussian regions we obtain below. We return to Eq. (|7]i 
when we examine the case of XENON IT in Sect. II V Cl 



IV. SIMULATION RESULTS 
A. The Maxwell-Boltzmann Case with Light WIMPs 

As a first example, we illustrate the importance of the un- 
certainty on the dispersion in setting exclusion limits. Null re- 
sults at MAJORANA Demonstrator, LUX and XENON IT 
would result in exclusion curves shown in Fig. [T] The solid 
curves assume the standard value vq — 230 km/s, while the 
bands show the result of varying vq within ±60 km/s. For 
comparison, the existing exclusion curves are also shown (for 
fixed Vq). As is clear from the Figure, the precise value of the 
velocity dispersion becomes particularly important for light 
WIMP masses, where the sensitivity of the experiments varies 
by an order of magnitude. Physically, this is due to the fact 
that the light WIMPs to which these experiments are sensitive 
all lie on the high-velocity tail of the distribution. The larger 



the velocity dispersion, the more stringent the bound for light 
masses. 

Let us now consider our main simulation, of an actual 
WIMP detection. We will examine in detail two useful bench- 
mark point shown in Fig. [T] one with input data coming from 
a 8 GeV WIMP with scattering cross section 10^^'^ cm^, 
and the other a 7 GeV WIMP with a 10^"'^ cm^ cross sec- 
tion. Both of these are just below the current XENONIOO 
bounds |21 1. At this mass and cross section, the MAJORANA 
Demonstrator would see nearly 500 events in the 8 GeV 
case and around 3300 events in 7 GeV case. 

For each of the two points, we generate "data", using — 
— 230 km/s and setting the escape speed to Vesc — 544 
km/s. We then perform a fit on a three-dimensional grid of 
points, in the (toxi ^^o) space. We keep v^g^ ~ 544 km/s 
and Ve = 230 km/s both fixed in the fit. This is done pri- 
marily for the sake of clarity: to avoid confusion among many 
parameters. 

Let us first examine the 8 GeV case. In Fig. |2]we show 
three slices of the full three-dimensional confidence regions at 
fixed cross sections. In addition to the input value (black), we 
also choose one lower value (in green) and one higher value 
(in red). Marginalizing (integrating likelihood) over the cross 
section yields the light blue regions. 

We see a clear degeneracy between mx and vq, which be- 
comes particularly prominent in the case of larger statistics (in 
the right panel). To understand it, let us specialize to the right 
panel in Fig [2] and consider a fixed cross section (the black 
region, for definiteness). Recall that for light WIMPS, v,nin is 
rather large (only high velocity WIMPS can produce enough 
recoil). This implies that the velocity integral g{vmin) (see 
the Appendix for the full expression) is well-approximated by 



g{Vmin) OC 



1 — erf 
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such that directions of constant {vmin—Ve)lvo are degenerate. 
To translate this into a curve in the mx — vq plane, we must 
fix Eb. or equivalently Vmin- Taking 3 keV falls roughly in 
the middle of the signal spectrum. 

We display the analytic estimate based on Eq. (|9]) as a 
dashed gray line in the right panel of Fig. [3] It nicely re- 
produces the shape of the black region in Fig. |2] Note that 
the contours at other cross sections can be understood with an 
analogous argument. 

Notice that the analytical contour bends away from the pre- 
diction once uo is greater than about 350 km/s. This is because 
of the effects of the finite the escape speed: for vq > 350 km 
and Vesc = 544 km/s , a significant part of the Maxwellian tail 
is cut off. (Eq. (j9]l is obtained neglecting this cutoff.) 

Let us show how this degeneracy impacts the ability of the 
experiment to determine the mass and cross section of dark 
matter To do so, we marginalize out the vq dependence by 
imposing a top hat prior on the velocity dispersion centered 
on the true value 230 km/s, for a variety of uncertainties ay 
(±25 km/s, ±75 km/s and from zero to the escape speed). The 
results for both examples are shown in the top two panels of 
Fig. [3] where we have also shown the regions coming from 
fixing the dispersion to 230 km/s (black). The right panel is 
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FIG. 2: Slices of the best-fit regions from the Majorana DEMONSTRATOR in {mx , <7n,vo) space. Light blue shows the cross section- 
marginalized regions, whereas the other contours are for fixed cross section as indicated in each plot. The dashed gray curve represent 
the analytic estimate, Eq. |9](, which agrees well with the black contour. The input data is generated from a Maxwellian halo with with 



vo 



vlsr ~ 230 km/s, with DM mass of 8 GeV and cross section 10 cm (left) and 7 GeV and cross section 10 cm (right). 



particularly instructive. We see that in the fixed astrophysics 
assumption one might expect an excellent measurement of the 
mass and cross section, but the measurement deteriorates sig- 
nificantly as the prior on get weaker and weaker. The long- 
tail extending to high masses results from the ability of high- 
mass/low-velocity dispersion cases to mimic the true data rea- 
sonably well. The low-mass tail, which originates from the 
ability of low-mass/high-dispersion cases to match the true 
data, is considerably shorter. 

Going one step further, we can marginalize out the cross 
section to see how well the WIMP mass, mx, can be deter- 
mined. The results are shown in the bottom panels in Fig. [2] 
We first examine the 8 GeV case. Interestingly, when the 
dispersion is varied, the experimental reconstruction prefers 
WIMP masses around 6 GeV. This is the result of projecting 
the "sock-like" region depicted in the panel above: the "an- 
kle" contains many points along the va direction, which when 
projected onto the mx axis, skews the likelihood toward low 
masses. 

The same basic degeneracy and the concomitant skewed 
likelihood functions carries over for the 7 GeV case, despite 
the large increase in statistics. In the right panel of Fig. [3] we 
repeat the above exercise for the example depicted in the right 
panel of Fig.|2] Again, this example includes a factor ^1 in- 
crease in signal events, which is still insufficient to completely 
remove the degeneracy. In fact, the relative effect of the de- 
generacy is especially significant as the 1 a region expands 
from mx = Ttp j GeV to GA^H GeV in going from fixed 
dispersion to complete uncertainty. 

In summary, by inferring the astrophysics (wq) from the 
data, we suffer a significant sensitivity loss to the light WIMP 



mass. 



B. Via Lactea II Velocity Distributions 

Do our findings depend on the particular choice of the mock 
spectrum? To check this, we repeat the analysis, but instead 
generate signal events from the high resolution N-body dark 
matter simulations borrowed from the Via Lactea II (VL-II) 
project ll58l . The effect of such distributions on direct detec- 
tion phenomenology has been previously studied in f59] l60l 
and constraints from existing experiments derived in [37l . The 
effect of tidal debris in VL-II was first observed in ||49l and 
subsequently applied to direct detection in [30^ '611, though 
here we focus on the full contents of the velocity distribution. 

As we will show, the existence of the mass-dispersion 
degeneracy does not depend on the input spectrum coming 
from a Maxwell-Boltzmann form of the velocity distribution. 
Moreover, an experimentalist is unlikely to be able to exclude 
the Maxwellian hypothesis if the true halo is of a Via Lactea 
form. Note that although the unsealed VL-II distribution is 
not a realistic approximation of the Milky Way halo, we in- 
clude it here for comparison. The scaled VL-II distribution is 
scaled such that its dispersion is 220 km/s. 

As can be seen in Fig.|4]for the 8 GeV case, the same qual- 
itative features of mass and dispersion mis-measurement per- 
sist as well, though differing somewhat in the quantitative de- 
tails. There we see the small effect induced from the raw VL- 
II distribution which favors low mass solutions compared to 
the other two distributions. An important observation is that 
the best-fit points have per degree of freedom that is very 
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FIG. 3: Marginalized 95% CL regions for the Majorana DEMONSTRATOR experiment with a 100 kg-yr exposure. In the upper panels we 
impose a top-hat prior on the dispersion, centered on the true value 230 km/s with varying degrees of uncertainty . In the bottom row, we 
project this into the one-dimensional mx distribution. Results from a (8 GeV, 10^*'' cm^) and (7 GeV, 10^^^ cm^) input spectra are depicted 
in the left and right columns respectively. 



close to one, indicating that the isothermal fit is good. In- 
terestingly, the values of mx and (T„ obtained in the best-fit 
point are ^15% off from the original inputs. Nevertheless, 
since in an actual experiment the true ("input") values are un- 
known, the fact that the distribution is non-Maxwellian will 
not be registered. 

Lastly, we ask how well direct detection can determine 
the dispersion of the halo. One way of addressing this is to 
marginalize over the mass and cross section to determine the 
dispersion likelihood function, shown in the bottom panel of 
Fig. [4] In the 8 GeV case (red curves) the dispersion like- 
lihood function from DEMONSTRATOR is extremely broad, 
with a poor ability to reconstruct the correct dispersion. Its 



best fit is peaked near 150 km/s, but with significant Ukelihood 
extending all the way to dispersions of 600 km/s. The high- 
statistics 7 GeV case fares a bit better with a best-fit value of 
173 km/s for the MB and a best-fit of 207 km/s for the scaled 
VL-11 distribution. Further, though there the 7 GeV starts to 
show visible separation of the three distributions, it is not suf- 
ficient for discriminating among velocity distributions. 



C. Heavier WIMPs and XENONIT 

In the above we have focused on the light WIMP case since 
the lowering of energy thresholds in upcoming experiments 
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FIG. 4: Here we compare the CL regions when for the same (mx = 
8 GeV, CT„ = 10"''^ cm^) example when coming from three differ- 
ent velocity distributions: Maxwell-Boltzmann (black, dashed), Via 
Lactea II raw distribution (red dotted) and Via Lactea II scaled to 
vo = 220 km/s (blue, solid). Bottom: Likelihood functions of the 
velocity dispersion for all three distributions from above for both the 
7 GeV and 8 GeV examples. Here solid, dashed, and dotted curves 
are MB, VL-II scaled and VL-II raw distributions. 



could easily lead to an immediate and huge signal. Significant 
statistics in the more canonical heavy WIMP case will require 
significantly larger experiments such as XENON IT i55il . 

We stress that the energy threshold we have chosen 
is conservative, and will hopefully be significantly lower 
given that XENON IT will have a larger light yield than 
XENON 100 155 1 . Moreover, if XENON IT (or LUX) can 
lower their energy thresholds, they will have unprecedented 
sensitivity to a hitherto unexplored region of WIMP parame- 
ter space. Like the XENONIO and XENONIOO experiments, 
XENON IT uses a combination of amplified charge signals 
(S2) and prompt scintillation light (SI) to discriminate signal 
from background. Such low-thresholds are for example poten- 
tially achievable with the use of an S2-only analysis such as 
the one already performed by XENONIO L20J . Alternatively, 



Xe, 2.2 ton-yr, 95% CL 




0.5 



100 200 300 400 500 
mx(GeV) 

FIG. 5: Here the signal spectrum comes from a vq = 230 km/s 
halo with WIMP mass 100 GeV and cross section 10"*^ cm^. In 
both figures, 68% (black) and 95% CL (red and blue) are displayed. 
Above: One slice of the best-fit region with the cross section fixed to 
the true value, but floating vo and mx - Bottom: The ensuing allowed 
regions coming from the canonical fit with fixed dispersion (black) 
and Vo marginalized contour. 



XENON IT may be able to use the ionized free electron sig- 
nal as their primary energy estimator in order to obtain energy 
threshold approaching a few keV ll62l . Using the outer lay- 
ers of the detector for self-shielding, the 3 tons of total liquid 
Xenon translate to 1.1 tons of fiducial detector mass. 

As an example, consider a 100 GeV WIMP with a lO"**^ 
cm^ scattering cross section, which produces ~'100 events at 



8 



our XENONIT mock-up.^ Though the relation describing the 
degeneracy at low masses (Eq. |9]) no longer holds, we see 
in the top panel of Fig. [5] that a similar degeneracy persists 
at XENONIT. That the analytic description of the degener- 
acy does not hold here should not be surprising since v,nin is 
rather small in this case. 

When the dispersion is marginalized out (bottom panel of 
Fig.|5]l, we see that while the dispersion has an impact on the 
parameter reconstruction, it is less significant than in the light 
WIMP case. In addition, the 95% CL contours (dashed) en- 
counter a well known high mass degeneracy |29|. This addi- 
tional degeneracy appears because at sufficiently high WIMP 
mass v„iin no longer contains any information about mx- 
Thus the only place in the differential scattering rate (Eq. |4]i 
containing the WIMP mass is in the ratio an /nix such that 
this combination forms an irreducible degenerate direction in 
parameter space. This example also offers poor prospects for 
dispersion determination, yielding at Ict, 144^^2^ km/s. 
This skewed parameter determination is again due to the high 
mass/low dispersion models seen in the top panel of Fig. |5] 
that well mimic the input spectrum. 

Lastly, one may wonder at what mass mx does the mass- 
degeneracy relation starts to take effect. In examining runs 
using our XENONIT setup, we find that by masses around 20 
GeV the skewed a — mx regions (such as Fig. [3]) and "sock- 
like" mx — vq regions (such as Fig. [2]l start to appear. Thus 
even for more canonical "high-threshold" experiments the de- 
termination of the dispersion can be critical for precise DM 
mass measurements. 



D. Streams 

Unbound streams of dark matter are generally expected 
in the hierarchical process of structure formation |63| as 
well as in non-standard halo models, such as the late-infall 
model |64|. The stream coming from the tidal disruption 
of the Sagittarius dwarf galaxy is perhaps the most well- 
studied such example Ii47ll48ll65ll66ll . though N-body simu- 
lations predict a wealth of velocity substructure. One of their 
most dramatic signatures is in annual modulation searches, in 
which they can disrupt the prediction of a sinusoidal variation 
in time Il4gll66l . 

However, even at the level of the nuclear recoil spectrum, 
dispersionless streams produce rather distinctive phenomenol- 
ogy, with the stream contribution to g{vmin) appearing as a 
relatively sharp step-function. Generally, one expects such 
streams to form a subdominant contribution to the local den- 
sity, ^ = Pst/ Pdm < 1. with the remainder of the local DM 
population described by an equilibrium distribution such as 
Maxwell-Boltzmann |67|. Thus suppose when faced with a 
signal of DM, one is interested in excluding the MB only hy- 



Note that here we use a "product of Poissons" likehhood function as in 
Eq. ^ since the number of signal events in each bin can be small, though 
our results to not deviate significantly from a analysis. 



pothesis. Then when fitting to a stream -i- MB signal with a 
MB distribution, we can roughly estimate an expected sen- 
sitivity to streams scaling as N^t > \/Nmb + ^bkg^ with 
NsIjNmb, and Nt,kg, being the number of events in the rele- 
vant energy bin(s) from the stream, MB halo, and background 
respectively. Thus, for irreducible backgrounds, the best sen- 
sitivity to streams occurs when the stream contributes to en- 
ergy bins where the background dominates the MB signal. 
There one expects the ability to exclude the MB-only hy- 
pothesis when Nst > \/Nbkg. We simulated spectra for the 
parameters of the Sagittarius stream, for which v^tr ~ 340 
km/s |48|. With this scahng, in the case of the 7 GeV, 10""*^ 
cm^ example studied eai'lier, the Majorana Demonstra- 
tor would allow 1(7 detection of the Sagittarius stream if its 
density is ^ > 6%. 

Let us consider the other extreme case, in which the back- 
ground is completely subdominant to the MB signal which has 
been recently studied in detail by {66\. There they have found 
that a 10 GeV WIMP capable of explaining the CoGeNT re- 
sults, would have a systematic mass underestimate and cross 
section overestimate 166|. The authors of |66|, find that a 5% 
Sagittarius stream could be detectable at the 2a level with a 
10 kg-yr exposure on the C-4 CoGeNT upgrade. We plan to 
return to this topic in the future, exploring the implications for 
the Majorana Demonstrator. 



V. DISCUSSION AND CONCLUSIONS 

At the outset of this paper we set out to ask a basic question 
that we may well be faced with shortly: If DM is a light WIMP 
with a large cross section, will we have sufficient information 
to infer both the particle physics and the local astrophysics of 
DM? We have found that an accurate reconstruction of both 
the mass and the dispersion will be unlikely even when the 
signal consists of many thousands of events. This is primar- 
ily due to a previously unexplored degeneracy between the 
WIMP mass and the local velocity dispersion. Indeed, the 
uncertainty on the average velocity of the local WIMP popu- 
lation is likely to be the largest factor inhibiting an accurate 
mass measurement. 

In addition, we have found that this conclusion holds if the 
signal spectra came from an isothermal halo or one described 
by the high resolution Via Lactea-II simulation. We stress 
as well that despite these simulations including tidal debris 
which is made up of unequilibrated structure 14911501 . an ex- 
periment like the MAJORANA DEMONSTRATOR will not be 
able to infer a deviation from the standard isothermal halo. 
The phenomenologically more distinctive streams, however, 
could be well-identified with sufficient statistics. Of course, 
improved sensitivity to streams and other tidal debris would 
come from the inclusion of modulation or directional infor- 
mation (for recent work see |68 -73 1), which we have ignored 
in this first analysis. 

In view of our findings, it is clear that dark matter mass de- 
termination could well hinge on our ability to better constrain 
the velocity distribution of the Milky Way by other means. 
Such an improvement could come about from improved ana- 
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lytic modeling, for example including self-consistent distribu- 
tions Ii32i . or from improvements in the numerical simulation 
of dark matter halos. We plan to return to these questions in 
future work. 
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Appendix: The Velocity Integral 

For completeness here we include a simple analytic expres- 
sion for the MB, f{v) cx e"" /"o velocity integral relevant for 



giVrmn) - ^^^^ 



where N = erf(xesc) - ^^p^e~^"<:. 
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